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ABSTRACT 

We present a 60 ks Chandra ACIS-S observation of the thermal composite 
supernova remnant 3C 391. The southeast-northwest elongated morphology is 
similar to that previously found in radio and X-ray studies. This observation 
unveils a highly clumpy structure of the remnant. Detailed spatially resolved 
spectral analysis for the small-scale features reveals that the interior gas is gen- 
erally of normal metal abundance and has approached or basically reached ion- 
ization equilibrium. The hydrogen column density increases from southeast to 
northwest. Three mechanisms, radiative rim, thermal conduction, and cloudlet 
evaporation, may all play roles in the X-ray appearance of 3C 391 as a "thermal 
composite" remnant, but there are difficulties with each of them in explaining 
some physical properties. Comparatively, the cloudlet evaporation model is fa- 
vored by the main characteristics such as the highly clumpy structure and the 
uniform temperature and density distribution over most of the remnant. The 
directly measured postshock temperature also implies a young age, ~ 4 x 10 3 yr, 
for the remnant. The postshock gas pressure derived from the NE and SW rims, 
which harbor maser spots, is consistent with the estimate for the maser regions. 
An unresolved X-ray source is observed on the northwest border and its spectrum 
is best fitted by a power-law. 



1. Introduction 

3C 391 (G31. 9+0.0), a supernova remnant (SNR) with irregular morphology, has been 
observed in several electromagnetic bands. Radio observations with the VLA reveal an 
elongated structure extending from the northwest (NW) to the southeast (SE), surrounded 
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by a shell except on the SE border, where it appears that the SNR has broken out of a dense 
region into an adjacent region of lower density (Reynolds & Moffett 1993). In X-rays, both 
Einstein (Wang & Seward 1984) and ROSAT (Rho & Petre 1996) observations show that 
the centroid of the soft X-ray emission sits in the SE region. The ASCA observation reveals 
a hard X-ray enhancement in the NW region, and confirms the decrease in hydrogen column 
density across the remnant from NW to SE (Chen & Slane 2001). Two hydroxyl radical 
1720 MHz maser spots are found along the radio shell (Frail et al. 1996). The masers and 
the strong enhancement in [01] 63/xm emission near the northwestern edge (Reach & Rho 
1996) both provide evidence for shock interactions with a molecular cloud. CO and other 
molecular line observations place the location of the remnant at the southwestern (SW) edge 
of a molecular cloud (Wilner, Reynolds, & Moffett 1998; Reach & Rho 1999). 

3C 391 is an example of the "thermal composite" or "mixed morphology" category of 
supernova remnants, which also includes W28, W44, IC 443, G349. 7+0.2, and others (Rho & 
Petre 1998). They generate bright thermal X-ray emission interior to their radio shells, and 
have faint X-ray rims. They are usually found to interact with adjacent molecular clouds, 
characterized by the hydroxyl radical maser emission (Green et al. 1997, Yusef-Zadeh et al. 
2003). Non-thermal diffuse X-ray emission inside SNRs is widely believed to be synchrotron 
radiation from pulsar wind nebulae, while the nature of the internal thermal X-ray emission 
seen in mixed morphology remnants is still uncertain. 

At least four distinct scenarios have been proposed to explain centrally-brightened X-ray 
morphology. The first scenario is radiative cooling of the rim gas. Under this hypothesis, the 
gas at the rim has been cooled down in the radiative stage of evolution, with a temperature 
so low (< 10 6 K) that its X-ray emission is very weak, while the gas in the inner volume is still 
hot enough to emit strong X-rays (e.g. Harrus et al. 1997, Rho & Petre 1998). The second 
mechanism invokes thermal conduction. It is suggested that thermal conduction in the 
remnant can smooth out the temperature gradient from the hot interior to the cooler shell, 
and increase the central density in response to the associated change in pressure conditions. 
This results in luminous X-rays in the interior, and is dominant in the radiative stage (Cox 
et al. 1999, Shelton et al. 1999). The third scenario invokes cloudlet evaporation in the SNR 
interior. When a SNR expands in an inhomogeneous interstellar medium (ISM) whose mass 
is mostly contained in small clouds, the clouds engulfed by the blast wave can be evaporated 
to slowly increase the density of the interior gas; as a result, the SNR appears internally 
X-ray brightened (White & Long 1991). The fourth suggestion is that the mixed morphology 
is a projection effect. For shell-like SNRs that evolve in a density gradient such as at the 
edge of a molecular cloud, if the line of sight is essentially aligned with the density gradient 
as well as the magnetic field, the SNRs will appear as thermal composites (Petruk 2001). 
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To date, X-ray studies have not yet provided strong evidence favoring any one of these 
mechanisms. Here we report on high resolution Chandra observations of 3C 391 aimed at 
addressing the nature of the centrally filled X-ray emission in this SNR. 

2. Observation and Data Analysis 

We observed SNR 3C 391 with the Advanced CCD Imaging Spectrometer (ACIS) on 
board the Chandra X-ray observatory on 03 August 2002 (Obs. ID 2786) for an exposure of 
62 ks. The target center (18 h 49 m 28 s .3, -00°56'15") was placed ~ 2' off the telescope aim- 
point for an optimal coverage of the supernova remnant on the S3 (#7) CCD chip (using 
very faint mode). 

The level 1 raw event data were reprocessed to generate a new level 2 event file to 
capitalize on the Chandra Interactive Analysis of Observations (CIAO) software package 
(version 2.3). In the reprocessing, bad grades have been filtered out and good time intervals 
have been reserved. We also removed time intervals with significant background flares (peak 
count rates ^ 3a and/or a factor of ^ 1.2 off the mean background level of the observation). 
This cleaning, together with a correction for the dead time of the observation, resulted a 
net 60.7 ks exposure for subsequent analysis. Furthermore, a correction was made for the 
time-dependence of the ACIS gain. 

2.1. Spatial Analysis 

For imaging analysis, we produced exposure maps in 0.3-1.5, 1.5-3.0, and 3.0-7.0 keV 
bands. These maps were used for flat-fielding, accounting for bad pixel removal as well as 
correcting for telescope vignetting and the quantum efficiency variation across the detector. 
The images in the three bands are adaptively smoothed (using CIAO tool csmooth with 
signal-to-noise ratio of 3) and exposure-corrected. The tri-color X-ray image (0.3-1.5 keV in 
red, 1.5-3.0 keV in green, and 3.0-7.0 keV in blue), overlaid with the 1.5 GHz radio contours 
(from Moffett & Reynolds 1994), is shown in Fig. la. Another tri-color image with 1.5 GHz 
emission in red, 0.3-3.0 keV in green, and 3.0-7.0 keV in blue is shown in Fig. lb. A broad 
band (0.3-7.0 keV) image overlaid with the radio contours are shown in Fig.lc. The two OH 
maser points (Frail et al. 1996) are labeled in these three maps. 

We searched for point-like sources in three broad bands, 0.3-1.5 keV (S), 1.5-7.0 keV (H), 
and 0.3-7.0 keV (B). A combination of source detection algorithms were applied: wavelet, 
sliding-box, and maximum likelihood centroid fitting (see Wang, Chaves, & Irwin 2003 for 
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details of the approach). The estimation of the count rate of a source was based on the data 
within the 90% energy-encircled radius (EER) determined with the calibrated point spread 
function of the instrument. The information obtained for the point-like sources detected this 
way on the S3 chip are summarized in Table 1. 

In order to produce a map of diffuse emission, we removed the point-like sources with 
each source region as a circle of twice the 90% EER. After filling the source-removed "holes" 
(in a swiss-cheese-like map) by interpolating the intensity from the surrounding regions, we 
produced a smoothed (with count-to- noise ratio of 6), exposure-corrected diffuse X-ray map 
in 0.3-7 keV (Fig. Id). Also shown in Fig. Id are the regions used for spectral fitting (see 
§2.2). 

These X-ray images display the SE-NW elongated morphology, similar to that observed 
in the earlier radio (Moffett, & Reynolds 1994) and X-ray (Rho & Petre 1996, Chen & 
Slane 2001) observations, but with considerable small-scale structure revealed by the superb 
angular resolution of Chandra. 

The X-ray maps reveal a highly clumpy structure of the remnant, with clumps or 
knots located in both the SE and NW parts. A bright hard enhancement is peaked at 
(18 h 49 m 27 s .O, -00°54'10") near the the NW border (see Fig. la); this corresponds to region 
#3 in Fig. Id. A complex mixture of knots is seen in the SE part of the remnant interior, in- 
cluding at least four bright enhancements in all colors centered at ~ (18 h 49 m 29 s .3, — 00°56'42 // ) 
(region #6), (18 h 49 m 32 s .9, -00°56'58") (region #7), (18 h 49 m 35 s .5, -00°57'04") (region #8), 
and (18 h 49 m 29 s .8, -00°57'45") (region #9). Similar clumpy structures are also seen in the 
interior of SNR N49B (Park et al. 2003). Several remarkable, very bright knotty features 
appear on the east and SE border of 3C 391. Also, the eastern knot seems to coincide with 
the location of the infrared source IRAS 18470-0100 (whose fluxes are 1.32, 4.66, 23.8, and 
236 Jy at 12, 25, 60, and 100 /zm, respectively). These bright knots on the border may be 
small clouds that have recently been shocked by the supernova blast wave. 

In the RGB 3-color map Fig. la, yellow and light-green colors indicate strong soft emis- 
sion in the SE part, while fainter soft emission appears brown in color and is seen in the 
middle section of the elongated body of the remnant. A faint arc appears in the rim of the 
northern brown patch. This arc/shell seems to connect with another shell like structure at 
the northeastern (NE) rim, which is evident in radio (Fig. lb) and discernible in X-rays, and 
harbors one of the previously detected OH masers (Frail et al. 1996). On the opposite side, 
no shell-like X-ray structure can be discerned along the southern and SW border, but faint 
diffuse X-ray emitting gas seems to spill outside the radio border in the southwest. 

In Fig. la, the green color representing the middle-band emitting gas almost pervades 
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the whole postshock region of the NW part. The pervasive green color in the NW, mixed 
with a small blue patch, shows hardening of the X-ray spectrum. This is consistent with 
gradual increase in hydrogen column density across the remnant (Chen & Slane 2001; see 
also §3.2 below). 

The radio contour overlays make it very apparent that there is a shock interaction with 
a molecular cloud on the NW border because of the flat morphology of the radio emission. At 
the west rim, a soft X-ray brightened slab appears to be very close to (just slightly behind) 
the radio peak emission (region #12), and may be related to a small dense region there. 

A few point sources are seen in the remnant's geometric center and the south and west 
portions (Fig.l and Table 1). The small number of counts collected make it difficult to 
determine their nature. A bright hard enhancement (18 h 49 m 22 s .3, — 00°53'34") appears as 
an unresolved source on the NW border (source #17 in Table 1; region #5 in Fig. Id). We 
discuss the spectrum of this source below (§3.3.). 

2.2. Spectral Analysis 

Taking advantage of Chandra's superb angular resolution, we have performed the spatially- 
resolved spectral analysis of 3C 391. Prior to the extraction of spectra, we removed the point- 
like sources detected on the S3 chip (see Table 1) except for weak sources with log(P) > —9.0 
(see Table 1 for the detailed definition of the false detection probability P) and those with 
a high possibility of being clumps of diffuse emission (judged by eye). We defined 13 small- 
scale regions, 3 medium-scale regions, and the region of the entire remnant (as diagramed in 
Fig. Id) for spectral investigation. Most of the small-scale regions are chosen to include the 
small features of X-ray enhancement such as the knots and the faint shell like structures. 
Two medium-size circles are used to compare the mean spectral variation between the inner 
SE and NW portions. The remaining area on the S3 chip outside region #17 for the entire 
remnant was used for background. For each chosen region, we use the CIAO script acisspec 
to extract the spectrum and obtain the weighted response matrix. The net count rates of 
these extracted spectra are tabulated in Table 2 and Table 3. Corrections for the continuous 
degradation in the ACIS quantum efficiency have been applied to auxiliary response files. 
All the spectra mentioned above were regrouped to contain at least 25 net counts per bin. 

There are three distinct line features, Mg Hea (~ 1.35 keV), Si Hea (~ 1.85 keV), and 
S Hea (~ 2.46 keV) in these spectra except the spectrum of a NW unresolved source (region 
#5), confirming their thermal origin (see Fig. 2 and Fig. 3). The Fe L complex could have 
significant contribution in the range ~ 0.5-1.5 keV. The emission diminishes rapidly above 
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the photon energy ~ 5 keV and no Fe Ka emission is observed. 

Using the Morrison & McCammon (1983) interstellar absorption, we test to fit the 
thermal spectra with various collisional equilibrium ionization and non-equilibrium ionization 
(NEI) models for the optically thin plasma of SNRs available in the XSPEC code. We find 
that our spectra can be best described by the VNEI model with NEI spectral version 2.0, 
which uses the Astrophysical Plasma Emission Database (APED) to calculate spectrum 
1 . Considering the evident line emission from Mg, Si, S, and possibly Fe L, we treat the 
abundances of these element species as free parameters. In most cases, however, thawing 
the abundances did not significantly improve the spectral fit relative to that using solar 
abundances. For a few cases, especially for the emission from medium-scale regions and from 
the entire remnant, thawing the abundances do indeed improve the fitting. The spectral fit 
results are tabulated in Table 2. The spectrum of the NW unresolved source (region #5) 
can be best fitted with an absorbed power-law emission (see the spectral fit parameters 
summarized in Table 3) and cannot be fitted with a blackbody emission. 

The spectral fits show that the diffuse emission from various regions have ionization 
parameters (n e t) close to or higher than 10 12 cm -3 s. This implies that the hot plasma in the 
SNR is very close to, or is basically in, the ionization equilibrium. 

The spectral fits also show that the diffuse emission from small-scale regions can be 
well fitted with solar abundances, and the emissions from larger regions can be fitted with 
abundances very close to solar values. 

The temperature of the gas interior to the SNR is~0.5-0.6 keV, with only small fluctua- 
tions. For instance, the lowest temperature (~ 0.46 keV) is found at the faint SE rim (region 
#13), and the highest temperature (~ 0.8 keV) is seen at the NE rim (region #10). 

The absorption column density is found to generally increase across the remnant from SE 
to NW. This is consistent with the conclusion of the ROSAT (Rho & Petre 1996) and ASCA 
studies (Slane & Chen 2001), and is also consistent with the emission hardening shown in 
the previous section. The absorption in the SE part is about ~ 2.9 x 10 22 cm -2 (region #14). 
The highest absorption (~ 4.1 x 10 22 cm" 2 s _1 ) is found in a NW region (region #4) which 
is in blue in the tri-color image (Fig. la). This hardest region was also seen in the ASCA 
study (see Fig. 3 in Chen & Slane 2001). A very high absorption (~ 3.7 x 10 22 cm _2 s _1 ) is 
also found at the westmost rim (region #12). The difference of mean N H between the central 
SE region (region #14) and the central NW region (region #15) is ~ 6 x 10 21 cm~ 2 . Such a 
variation in N H is consistent with the existence of a molecular cloud to the NW (Wilner et 
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al. 1998; Reach & Rho 1999). 

3. Discussion 

Similar to the hydrogen column density distribution across the remnant obtained in 
previous X-ray studies, the N H values obtained in the spectral analysis are confirmed to 
increase from the SE toward the NW. The spatial analysis also shows an obvious spectral 
hardening in the NW part. This could naturally be explained by the scenario that the 
NW part of the remnant is embedded in the molecular cloud (Reynolds & Moffett 1993). 
The mean difference of N H between the SE and NW interior portions, ~ 6 x 10 21 cm -2 , is 
consistent with that found in the ASCA study, and implies a mean density of the ambient 
molecular gas (n(B.2)) of order ~ 10-20 cm -3 on the assumption that the SNR line-of-sight 
depth in the cloud is similar to the remnant radius (Chen & Slane 2001). We also adopt the 
distance d ~ 8d s kpc estimated there. 

According to the NEI fit to the spectrum of the entire remnant (region #17), the X-ray 
luminosity in 0.5-10 keV of the remnant is Lx = 3.5 x 10 36 <ig ergss -1 . From the volume 
emission measure frien^V ~ 1.5 x 10 59 rfg cm -3 , where / is the filling factor of the X-ray emit- 
ting gas, the mean number density of the gas interior to the remnant is 1.9/ -1 / 2 dg cm -3 
(here n e ~ 1.2n H has been used and a cylinder-like volume of a diameter 5' and a height 7' 
has been assumed). The remnant volume is V ~ 5.1 x 10 58 rf| cm" 3 . The mass of the X-ray 
emitting gas is M x ~ lAn H m H fV ~ 114/ 1 / 2 rfg /2 M . 

Raw estimates for the gas density of the regions used for spectral fit (apart from #17, 
the entire remnant) can also be given in terms of the volume emission measures. We assume 
that the circular regions correspond to spheres, the elliptical regions to oblate spheroids, and 
the rectangular regions to shell structures. The gas densities thus obtained are listed in the 
last column of Table 2. The densest X-ray emitting gas is located in the SE and eastern 
regions #1 and #2 that contain very bright knots (with density ^ 10/ -1 / 2 <i 8 cm -3 ). The 
other bright knots (regions #3 in the NW and #6, #7, #8, and #9 in the inner central SE 
portion) have a gas density ~ 5-7/ -1 / 2 rf 8 1//2 cm -3 . The regions along the remnant border 
(#10, #11, #12, #13, and #16) have a density ~ ISf-^dg 1/2 cm -3 . 

3.1. Origin of centrally-filled X-ray morphology and dynamics of the remnant 

As mentioned in the Introduction, four main mechanisms have been proposed to explain 
the X-ray morphology of thermal composite SNRs: projection effect, radiative shells, thermal 
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conduction, and cloudlet evaporation. Here we compare them with the properties found from 
our spatially-resolved spectral analysis. 

3.1.1. Projection effect 

The projection effect does not match the properties found in 3C 391. As mentioned 
above, the molecular cloud is located in the N W but the X-ray emission is not only enhanced 
in the NW half, but also in the SE half. Moreover, the hydrogen column density increases 
along the SE-NW elongation direction, so the density gradient of the ambient medium seems 
to be close to the projection plane, and hence will make the projection effect weak according 
to the Petruk (2001) model. 

3.1.2. Radiative rim 

Though there is no report in the available literature of the presence of an HI shell around 
the remnant like that in W44, another thermal composite SNR (Koo & Heiles 1995), the 
diffuse filamentary near-infrared [Fe II] and the mid-infrared 12-18 fim [Ne II] and [Ne III] 
emission along the NW radio shell provide some evidence for a radiative cooling at the rim 
(Reach, Rho, & Jarrett 2001). As viewed in Fig.l, the volume of the remnant as delineated by 
the radio shell (including the regions of the near- and mid-infrared emission) is basically filled 
with X-ray emitting gas. Several examined portions of X-ray emission on the borders, such 
as the NE rim (region #10), the northern rim (region #11), the western rim (region #12), 
the SW faint rim (region #16), and the SE rim (region #13), which should be dominated 
by the postshock gas of the supernova blast wave (and are indeed located at the edge of the 
radio shell), are found to be at temperatures higher than 5 x 10 6 K. This indicates that a 
considerable amount of gas at the blast shock has not yet suffered heavy radiative cooling, 
and is capable of generating considerable X-ray emission along the remnant's border. If the 
gas density of these regions, ~ l-3/ -1 / 2 G?g cm -3 , is approximated as the postshock density, 
then the preshock intercloud medium (ICM) has a density about 0.2-0.7/ _1 / 2 rf 8 ^ 2 cm~ 3 . 

According to the infrared spectroscopic study, the ambient environment in the NW part 
has complicated physical properties, with at least three typical preshock components (Reach 
& Rho 2000). For the atomic component implied by the [OIII] and [NIII] emission, the 
preshock density n a < 1 cm -3 and the shock velocity v s ~ 500 km s -1 ; for the moderate- 
density molecular component revealed by the [OI], [Sill], and [Fell] emission, the gas density 
n c i ~ 10 2 cm~ 3 and the shock velocity v d ~ 10 2 kms _1 ; and for the dense clumps giving 
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rise to H 2 , OH, CO, CS, and H 2 molecular emission, the gas density ~ 10 4 cm" 3 and the 
shock velocity ~ 20kms _1 . Thus there may be two possibilities about the coexistence of the 
radiative filaments and the X-ray emitting gas along the NW rim: 1) the cooling shell may 
be in the early stage of its formation; or 2) the radiative filaments may be caused by the 
propagation of the shock waves into some regions of higher density in the molecular cloud 
so that these regions of postshock gas have cooled down severely while most of the gas at 
the rim is still hot. 

In the first case, we note that the bright [Fell] emission coincides precisely with the 
brightest radio bar along the western border, which is regarded as the interface between the 
remnant and the molecular cloud (Reach, Rho, & Jarrett 2001), and that, as we find here, 
the soft X-ray brightened slab along the west border (in region #12) appears to be very close 
to, but slightly behind, the radio bar. The shell formation time can be estimated according 
to tshdi « 5.3 x 10 4 (^ 14 /no /7 ) yr (Cox et al. 1999). If the mean preshock hydrogen nucleus 
density is n ~ 30 cm" 3 (twice the molecular density (n(H 2 )) ~ 10-20 cm" 3 given above), 
^sheii ~ 7.6 x 10 3 i?5( 14 yr; if the preshock density is n ~ 100 cm" 3 for the moderately dense 
cloud, t s heii ~ 3.8 x 10 3 i?5i 14 yr. The remnant's age would not be much larger than these 
estimates. 

In the second case, the blast shock propagates in an inhomogeneous medium. After 
a blast shock impacts a cloud, a transmitted shock (at velocity v c i) is expected to move 
into the cloud. The pressure balance between the shocked cloud and the shocked ICM gives 
n ci v d = ^ n icmV 2 s i where factor (3 ~ 1 for a shock interaction with small clouds (McKee & 
Cowie 1975) and j3 ~ 4.4 and 6 for a strong shock hitting a rigid cloud plane for n c i/n = 100 
and oo, respectively (Zel'dovich & Raizer 1967). The moderate-density (n c / ~ 100 cm" 3 ) 
cloud may, after it is shocked (with v d ~ 100 cm" 3 ), be responsible for the filamentary 
near- and mid-infrared emission and the tenuous atomic component {n icm < 1 cm" 3 ) may, 
after shocked, correspond to the postshock X-ray emitting gas detected by Chandra. For 
the western X-ray slab (kT ~ 0.6 keV) along the radiative filaments, the velocity of the 
blast wave in the hot gas is v* = (16kT/3fim H ) 1/2 ~ 700 kms" 1 (where the mean atomic 
weight \x = 0.61), and the preshock density is nj cm = n H /4 ~ 0.5 cm -3 . Thus we have 
(3 ~ 4.1/ 1 / 2 d<!/ 2 , which suggests that the SNR shock wave hits the plane surface of a dense 
cloud. This is consistent with the flat morphology of the NW radio bar. 

3.1.3. Thermal conduction 

The spatially- resolved spectral analysis demonstrates that the hot gas over most interior 
regions of the SNR is at similar temperatures of around 0.5-0.6 keV. Such a uniform distri- 
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bution of the gas temperature inside the remnant is not consistent with the canonical Sedov 
distribution in which the inner temperature is higher than that near the rim. On the other 
hand, the X-ray morphology of this remnant is not limb brightened, and this is inconsistent 
with the Sedov solution as well. Cox et al. (1999) explain this sort of morphology (as ob- 
served in W44) in terms of the effect of thermal conduction of the interior gas, which could 
prevent the formation of a hot vacuous cavity described by the adiabatic Sedov solution and 
provide enough material in the center to make it bright in X-rays. 

The thermal conduction scenario could apply to 3C 391 to some extent. If the sup- 
pression of conduction by magnetic fields can be ignored, the conduction timescale would be 
£ C ond ~ n e k£ 2 / k, where £ denotes the linear scale of the temperature gradient and the conduc- 
tivity is given by k — 1.84 x 10 _5 T 5//2 /lnAergss _1 K _1 cm -1 with the Coulomb logarithm 
In A = 29.7 + In n7 1/2 (T/10 6 K) (Spitzer 1962). If £ is taken to be of order the remnant's 
mean radius 7pc, the timescale would be t con d ~ 4.7 x 10 4 yr, much larger than the rem- 
nant's inferred age. It would, however, be much smaller than this in the early stage of the 
evolution in which T would be higher and £ would be smaller. If we look at small spatial 
scales of order the separation between clumps (typified by 1' ~ 2.3 pc), the timescale would 
be 5.2 x 10 3 yr, comparable to the remnant's age (see §3.1.4.), and hence thermal conduction 
may play a role in smoothing of the interior temperature profile. We note that X-ray clumps 
are also present in W44 (Shelton, Kuntz, & Petre 2004) where thermal conduction may be 
an efficient mechanism. The fact that the hot gas in 3C 391 is near ionization equilibrium 
suggests that there may be relatively little newly shocked material and is consistent with the 
idea of some mechanism smoothing the properties. More significantly, the presence of the 
radiative filaments along the NW border is also in favor of the thermal conduction model. 

There are, however, disagreements in the distribution of physical properties between the 
observation and the model interpretation. In the thermal conduction scenario in which a 
centrally brightened X-ray morphology is reproduced (Cox et al. 1999, Shelton et al. 1999), 
the central pressure is about 0.3 of that at the edge and the central density would be about 
much lower (0.13 times at radiative shell formation) than the preshock ICM density, but this 
is not consistent with the case of 3C 391. In fact, this remnant has a uniform temperature 
around 0.5-0.6 keV throughout and the inner mean density ~ 3f- l/2 d% 1/2 cm 3 (regions #14 
and #15) is not smaller than the hot gas density ~ l-3/ -1 / 2 d g cm' 3 along the periphery. 
In the model the temperature would be expected to decrease gradually away from the inner 
portions to the outer. Unlike W44 in which the gas temperature is found to drop by a factor 
^2 between the center and ~ 6pc from the center (Shelton et al. 2004), in 3C 391 our 
measurement does not suggest the temperature decrease even if the line-of-sight projection 
is taken into account. The inner temperature is not higher than the gas temperature in most 
portions along the rim, and highest temperatures appear at a location on the border (region 
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#10, the NE rim). 

3.I.4. Cloudlet evaporation 

The phenomenon of low temperature in the inner portions could result from substantial 
cloudlet evaporation because part of the thermal energy is deposited to the gas evaporated 
from the cold cloudlets. The relatively uniform distribution of temperature (even with 
slightly lower values at the center) is actually expected by the White & Long (1991) cloudlet 
evaporation model for model parameters r — > 00 and C/r 3 (see Fig. 2 and Fig. 4 therein), 
where r is the ratio of the cloud evaporation timescale to the SNR's age and C is the 
ratio of the mass in the cloudlets to the mass of ICM. The ratio between the mean density 
(~ 2f~ 1 / 2 d 8 1 ^ 2 cm~ 3 ) and the density along the border (~ l-3/ _1 / 2 <i 8 1 ^ 2 cm~ 3 ) is basically 
consistent to that predicted in the evaporation model. In fact, the highly clumpy structure 
unveiled in this observation lends support to the conjecture that the ambient molecular cloud 
is inhomogeneous. The cloudlets engulfed by the supernova blast wave can act as a large 
reservoir of interior gas by gradual evaporation. This could also explains why the mean 
density of the hot gas inside the remnant (as mentioned above) is much lower than that of 
the ambient cloud gas. The association between the thermal composites and the OH maser 
emission in molecular clouds is strongly suggestive of the role played by clouds in the SNRs 
(Yusef-Zadeh et al. 2003). 

With a little adjustment of the model parameters used in the ASCA study of 3C 391 
(Chen & Slane 2001, see Table 6 therein), we here adopt C/r ~3-6. Using the temperature 
measured from the entire remnant (region #17) kTx ~ 0.56 keV as the average tempera- 
ture, the White & Long (1991) model would predict a postshock temperature 0.46-0.64 keV, 
consistent with most of the postshock temperatures measured from the spectral analysis. 
The velocity of the blast wave is v s = (l6kT s /3u,m H y/ 2 ~ 620-730 kms" 1 (where the mean 
atomic weight \i = 0.61). Adopting a mean radius of the remnant r s = 3'r 3 ~ 7.0r 3 d 8 , 
we estimate the dynamical age of the remnant: t = 2r s /5v s ~ 3.7-4.4 x 10 3 r 3 (i 8 yr. Be- 
cause our Chandra observation can measure the postshock temperature of 3C 391, this age 
estimate should be more accurate than those obtained from ROSAT and ASCA observa- 
tions using the same model, and the age of remnant is somewhat smaller than previously 
thought. The preshock gas density (n ~ 0.2-0.7/ _1 / 2 dg 1//2 cm" 3 ) estimated above could be 
regarded as the preshock ICM density. An alternative estimate of the preshock ICM density 
is obtained from the remnant's X-ray luminosity [Lx = 3.5 x 10 36 <ig ergss -1 ): n ~ 0.07- 
0.4 cm -3 by use of eq.(21) in White & Long's (1991) model. According to this model, the 
explosion energy of the supernova remnant is E = [167r(1.4n m H )/25(7+ l)K](rl/t 2 ) ~ 0.3- 
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1.4 x 10 51 r|<ig(n /0.3 cm 3 ) ergs (where the adiabatic index 7 = 5/3 and the ratio of thermal 
to kinetic energy K ~ 0.385-0.132). 

The application of the evaporating cloudlets model to 3C 391 may raise a few questions, 
which are discussed below and need further clarification. 

First, the remnant ages may be systematically underestimated by the White & Long 
model for a given explosion energy. W44 is a typical example, where an age of ~ 7 x 10 3 yr 
estimated from the evaporation model (Rho et al. 1994) is smaller than the spindown time 
scale ~2x 10 4 yr of the associated pulsar (Wolszczan, Cordes, & Dewey 1991). The age 
derived above for 3C 391 (~4x 10 3 r 3 rf 8 yr) using the evaporation model may possibly be an 
underestimate. If a Sedov evolution is adopted, the postshock temperature kT s ~ 0.6 keV 
would correspond to a similar age 2r s /5v s ~ 4 x 10 3 r 3 dsjr; as a comparison, Reynolds & 
Moffett (1993) present an estimate ~ 1.7 x lO 4 ^^ 1 yr for a cloud density n ~ 100 cm -3 . 
Considering the possibility of the SNR being in the radiative phase, Chen & Slane (2001) 
obtain an age ~ 1.9 x lO^E^ 42 (r^dg) 10 / 3 yr for a preshock density n ~ 30cm~ 3 . 

Second, the cloudlet evaporation scenario needs a model parameter r — > oo (with 
C/t^3), and thus the influence of cloud evaporation on the properties of the interior gas 
seems possibly negligible. Similar parameters are required to explain the X-ray morphology 
of another thermal composite MSH11-61A observed with ASCA (Slane et al. 2002). In con- 
trast, hydrodynamic simulations (e.g., Xu & Stone 1995, using a density contrast x — 10) 
show that clouds tend to get disrupted soon after passing through a strong shock. In our 
observation of 3C 391, however, the SE bright knots (in region #1) appear compact and sur- 
vive the shock passage. The ratio of the hot gas density of these knots to the mean density 
of the entire remnant is about 10 (see Table 2), so it would be reasonable to assume that the 
cloud-to-ICM density contrast x is much higher than this number or to assume a high mass 
contrast C. The other X-ray knots interior to the remnant may also be the dense regions 
that have not been destroyed by the blast shock. The denser clouds could be expected to 
survive longer time, and still a large amount of dense matter expected to have not been evap- 
orated from cloudlets yet. Infrared observations of 3C 391 have detected dense ( ^ 10 4 cm~ 3 ) 
molecular clumps in the southwest "broad molecular line" region, each H 2 clump being of 
size ~ 0.1 pc (Reach & Rho 1999, Reach & Rho 2000, Reach et al. 2002). (As a comparison, 
the compact SE knots in region #1 are of size ~ 0.5 pc.) It is not impossible that the X-ray 
knots contain similar dense cloudlets. 

Third, if the bright knots stands for the newly evaporated cloud material, it should 
have been seen cooler than the mean temperature. Nonetheless, because the evaporated gas 
would be soon heated up to the environmental temperature about 2 times the cloud radius 
from the cloud (Cowie & McKee 1977), the spectral properties of the cooler component could 
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be smoothed out if the knots are much larger than the cold cloudlets (for instance, 0.5 pc 
versus 0.1 pc as described above). 

Fourth, the newly evaporated gas should have a low ionization age. The knots in 
regions #2, #6, #7, #8, and #9 marginally satisfy this condition (>10 n s), especially if 
the evaporated gas had been pre-ionized by the ultraviolet photons and/or the transmitted 
shock wave, while those in region #1 has a high ionization age (> 10 12 s). The latter case 
does not satisfy the condition, unless, again, the low ionization age could be smeared by the 
surrounding gas if the cold clouds are very small. 

Concluding this section, we suggest that three mechanisms, radiative rim, thermal con- 
duction, and cloudlet evaporation, all play roles in the inner brightened X-ray morphology 
of 3C 391. By comparison, the cloudlet evaporation can largely explain some important 
properties such as the highly clumpy structure, the uniform temperature and density distri- 
bution almost all over the remnant, and the much lower X-ray emitting gas density than the 
mean ambient cloud density. Nonetheless, even cloudlet evaporation is driven by saturated 
conduction. It is probably true in many thermal composite remnants, and in 3C 391 as well, 
that the morphology arises from a combination of expansion into an inhomogeneous medium 
and heat conduction on various scales. 



3.2. Maser spots 

As in many other thermal composites, hydroxyl radical 1720 MHz maser spots are found 
within 3C 391. The two maser spots are located on the radio shell (Frail et al. 1996), one 
of which is here found to be harbored in the NE rim (region #10) where X-ray emission is 
modestly enhanced at a high plasma temperature (~ 0.8 keV) and the other of which is in 
the southwest region (#16) where the X-ray emission is faint (see Fig.l). 

The OH (1720 MHz) maser spots are believed to be due to C-type shock collisional 
excitation in the clumpy molecular clouds, with cloud temperature and density conditions 
50 < T < 125 K and 10 4 < n H2 < 5 x 10 5 cm" 3 (Elitzur 1976; Lockett, Gauthier, & Elitzur 
1999). The production of OH molecules is enhanced by the X-ray emission (with an ionization 
rate ^lO^s" 1 ) behind the C-shock (Wardle 1999). 

According to Wardle's model formula, the X-ray induced ionization rate in the post- 
shock regions of 3C 391 is about 2 x 10~ 15 s" 1 (using the X-ray luminosity Lx ~ 2.3 x 
10 36 ergss _1 and the mean radius ~ 7pc), well above the rate needed for the OH enhance- 
ment. Let us take the X-ray enhanced region #10 and X-ray faint region #16 as oblate 
spheroids. Thus the volume in region #10 is ~ 6.3<i|pc 3 and the mean gas density is 
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~ 2.3f~ 1 / 2 d 8 1 ^ 2 cm -3 . So, the pressure in the hot intercloud gas, p icm 2.3n H kT ~ 6.7 x 
lO" 9 /" 1 / 2 ^ 1//2 ergs cm" 3 . Similarly, for region #16, the volume is ~ 240rf|pc 3 , the mean gas 
density is ~ l.l/ -1 / 2 d 8 ^ 2 cm" 3 , and the gas pressure is p icm ~ 2.1 x lO" 9 /" 1 / 2 ^ ^ 2 ergs cm" 3 . 
Such gas pressures are consistent with the that estimated for the maser portions, p c \ = 
n m kT ~ 0.07-8.6 x 10" 9 ergs cm" 3 . The far-infrared H2O and OH emission lines are found 
around the SW maser site. They are consistent with the passage of shock wave through dense 
clumps and the postshock gas is estimated to have a density ~ 2 x 10 5 cm" 3 and temperature 
100-1000 K (Reach & Rho 1998). Such density and temperature are similar to the condi- 
tion for the maser production, with a pressure range p c \ ~ 2.8 x 10~ 9 -2.8 x 10~ 8 ergs cm" 3 , 
comparable to the above pressure estimates for the SW maser region. 

3.3. Stellar remnant? 

The association of 3C 391 with a dense molecular cloud makes it possible that 3C 391 is 
the interstellar remnant of the supernova explosion of a massive star, which may also leave 
behind a compact star as a result of gravitational core collapse. It is thus important to see 
if there is any evidence for this stellar remnant. We do locate numerous point-like sources 
within the boundary of the remnant (for example, evident point sources #22 and #26 right 
in the center, #24 and #32 in the south, etc.) (see Fig. Id and Table 1). However, spectral 
analysis for these sources cannot be carried out because of the small number of counts. 

A bright unresolved source (source #17) appears near the NW rim and seems to be 
located at the north apex of the NW diffuse hard emission (the blue patch in Fig. la and 
Fig. lb). The spectrum of this source extracted from region #5 (with 158 net counts) is 
better described with a powerlaw than with a bremsstrahlung or a blackbody (see Fig. 4 and 
Table 3). The hydrogen column density obtained is similar to those of other portions of 
this remnant, so the possibility of the association of this source with the remnant cannot 
be ruled out. If this source is at the distance of 3C 391, its unabsorbed X-ray luminosity 
is ~ 7.0 x lO^dsergss -1 . The steep photon index of the possible powerlaw model 2.5-4.5 
is higher than the typical indices 1.2-2.2 of active galactic nuclei (Turner & Pounds 1989), 
but is similar to those of the some "compact central objects" within SNRs (e.g. Pavlov et 
al. 2002). If the association with the gaseous remnant is true, the stellar remnant may have 
moved away from its explosion site. If we follow Reynolds & Moffett's assumption that the 
explosion site is located at the center of the NW half, then the source would have displaced 
2' toward the border. Taking an age of ~ 4 x 10 3 yr, then the displacement velocity would 
be ~ 1100 km s" 1 . Such a transverse velocity is quite high, but not unreasonable; pulsar 
proper motion studies have provided evidence that a small population of neutron stars have 
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velocities in excess of ^ 1000 km s" 1 (Lai 2003). For example, the pulsar B2224+65 in the 
Guitar nebular pulsar has a velocity ~ 800-1600 kms -1 (Chatterjee & Cordes 2004). Since 
the age used may be an underestimate, a higher remnant age would allow for a lower velocity. 

If the unresolved source is an offset pulsar, a trailing pulsar wind nebula might be 
expected to be present to the south, roughly in the hard emission region (the blue patch 
mentioned in §2.1). However, the spectrum of the diffuse gas of the blue patch (region #4) 
is well described by a thermal model (Table 2). If a power law component is added to fit 
the spectrum, the photon index would be 7.71 2 .' \, too steep to be physically true. Since the 
spectrum of region $4 is very week below 1.3 keV (Fig. 2) and the hydrogen column density 
is very high (~ 4.1 x 10 22 cm -2 ), this hard region seems to be the diffuse hot gas suffering 
heavy absorption in soft X-rays. 

4. Conclusion 

We observed the thermal composite supernova remnant 3C 391 for 60 ks using Chandra 
ACIS-S detector and carried out a detailed spectro-spatial X-ray analysis which results in 
following conclusions. 

1. The southeast-northwest elongated morphology is similar to that previously found in 
radio and X-ray studies. Faint shell-like X-ray structures appear along the northern 
and NE boundary and a faint slab appears at the west rim, while on the SW side the 
diffuse gas seems to expand out of the radio boundary. The X-ray emission is hardened 
in the NW, consistent with the increase of hydrogen column density from SE to NW, 
and all the images and spectra show compellingly that the intervening column density 
is the primary influence on the broad band X-ray appearance of the remnant. 

2. This observation unveils a highly clumpy structure in the remnant. The spatially re- 
solved spectral analysis for the small-scale features found that the gas is generally of 
normal metal abundance and has approached or basically reached ionization equilib- 
rium. 

3. Three mechanisms, radiative rim, thermal conduction, and cloudlet evaporation, may 
all play roles in the X-ray visage of 3C 391 as a "thermal composite" SNR, but there are 
difficulties with each of them in explaining some physical properties. By comparison, 
the cloudlet evaporation model can largely explain some important properties such as 
the highly clumpy structure, the uniform temperature and density distribution almost 
all over the remnant, and the much lower X-ray emitting gas density than the mean 
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ambient cloud density. The enhanced emission in both the SE and NW halves as well 
as the hydrogen column density gradient does not support the projection effect model. 
The directly measured postshock temperature implies a young age, ~ 4 x 10 3 yr, for 
the remnant. 

4. The postshock gas pressure derived from the NE and SW rims, which harbor maser 
spots, is consistent with the estimate for the maser regions. 

5. An unresolved source on the northwest border is best fitted in spectrum by a power- 
law, and its position and absorption column density might suggest a possibility for an 
association with 3C 391. 



YC and YS acknowledge support from NSFC grants 10073003 & 10221001 and grant 
NKBRSF-G19990754 of China Ministry of Science and Technology, POS acknowledges sup- 
port from NASA contract NAS8-39073 and grant GO2-3081X, and QDW acknowledges 
NASA LTSA grant NAG5-8935. The authors thank an anonymous referee for helpful com- 
ments and advices. 



-17- 



REFERENCES 

Chatterjee, S. & Cordes, J.M. 2004, ApJ, 600, L51 
Chen, Y. & Slane, P.O. 2001, ApJ, 563, 202 
Cowic, L.L. & McKee, C.F. 1977, ApJ, 211, 135 

Cox, D. P., Shelton, R. L., Maciejewski, W., Smith, R. K., Plewa, T., Pawl, A., & Rozyczka, 
M. 1999, ApJ, 524, 179 

Elitzur, M. 1976, ApJ, 203, 124 

Frail, D.A., Goss, W.M., Reynoso, E.M., Giacani, E.B., Green, A.J., & Otrupcek, R. 1996, 
AJ, 111, 1651 Goss, W.M., Winkler, P.F., & Williams, B.F. 1997, AJ, 113, 1379 

Green, A. J., Frail, D. A., Goss, W. M., Otrupcek, R. 1997, AJ, 114, 2058 

Harms, I. M., Hughes, J. P., Singh, K. P., Koyama, K., & Asaoka, I. 1997, ApJ, 488, 781 

Koo, B.-C. & Heiles, C. 1995, ApJ, 442, 679 

Lai, D. 2003, in 3D Signatures of Stellar Explosion, a workshop honoring J.C. Wheeler's 
60th Birthday (astro-ph/0312542) 

Lockett, P., Gauthier, E., & Elitzur, M., 1999, ApJ, 511, 235 

McKee, C.F. & Cowie, L.L. 1975, ApJ, 195, 715 

Moffett, D. A., & Reynolds, S. P., 1994, ApJ, 425, 668 

Morrison, R., & McCammon, D., 1983, ApJ, 270, 119 

Park, S., Hughes, J. P., Slane, P.O., Burrows, D.N., Warren, J.S., Garmire, G.P., & Nousek, 
J.A. 2003, ApJ, 592, L41 

Pavlov, G.G., Sanwal, D., Kiziltan, B., & Garmire, CP 2001, ApJ, 559, L131 

Petruk, O. 2001, A&A, 371, 267 

Reach, W.T., & Rho, J.H., 1996, A&A, 315, L277 

-, 1998, ApJ, 507, L93 

-, 1999, ApJ, 511, 836 



-, 2000, ApJ, 544, 858 
Reach, W.T., Rho, J.H., & Jarrett, T.H., 2002, ApJ, 564, 302 
Reynolds, S.R, & Moffett, D.A. 1993, AJ, 105, 2226 
Rho, J.H., & Petre, R. 1996, ApJ, 467, 698 
Rho, J.H., & Petre, R. 1998, ApJ, 503, L167 

Rho, J.H., Petre, R., Schlegel, E.M., & Hester, J.J. 1994, ApJ, 430, 757 

Shelton, R. L., Cox, D. P., Maciejewski, W., Smith, R. K., Plewa, T., Pawl, A., & Rozyczka, 
M. 1999, ApJ, 524, 192 

Shelton, R. L., Kuntz, K. D., & Petre, R. 2004, ApJ, accepted, astro-ph/0407026 

Slane, P., Smith, R., Hughes, J. P., & Petre, R. 2002, ApJ, 564, 284 

Spitzer, L. Jr. 1962, Physics of Fully Ionized Gases (New York: Interscience) 

Turner, T.J. & Pounds, K.A. 1989, MNRAS, 240, 833 

Wang, Q.D., Chaves, T., & Irwin, J. A. 2003, ApJ, 598, 969 

Wang, Z.R. & Seward, F.D. 1984, ApJ, 279, 705 

Wardle, M. 1999, ApJ, 525, L101 

Wilner, D.J., Reynolds, S.P, & Moffett, D.A. 1998, ApJ, 115, 247 
White, R.L., & Long, K.S. 1991, ApJ, 373, 543 
Wolszczan, A., Cordes, J.M., & Dewey, R.J. 1991, ApJ, 372, L99 
Xu, J. & Stone, J.M. 1995, ApJ, 454, 172 

Yusef-Zadeh, F., Wardle, M., Rho, J., & Sakano, M. 2003, ApJ, 585, 319 

Zel'dovich, Ya. B. & Raizer, Yu. M. 1967, Physics of shock waves and high-temperature 
hydrodynamic phenomena (New York: Academic Press) 



This preprint was prepared with the AAS IATgX macros v5.2. 



-19- 



Table I. Chandra List of Point-like Sources on S3 Chip 
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Table 1 — Continued 
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Note. — Column (1): Generic source number. (2): Chandra X-ray Ob- 
servatory (unregistered) source name, following the Chandra naming conven- 
tion and the IAU Recommendation for Nomenclature (e.g., http://cdsweb.u- 
strasbg.fr/iau-spec.html). (3): Position uncertainty (la) in units of arcsec. 
(4): The false detection probability P that the detected number of counts 
may result from the Poisson fluctuation of the local background within the 
detection aperture [log(P) smaller than -20.0 is set to -20.0]. (5): On-axis 
(exposure-corrected) source count rate in the 0.3-7 keV band. (6): The hard- 
ness ratio defined as HR = (H — S)/(H ± S), where S and H are the net source 
count rates in the 0.3-1.5 and 1.5-7 keV bands, respectively. The hardness 
ratios are calculated only for sources with individual signal-to-noise ratios 
greater than 4 in the broad band (B=S+H), and only the values with uncer- 
tainties less than 0.2 are included. 

a These sources are not removed for spectral analysis, either because of 
high log(P) (we define a threshold > —9.0) or high possibility of being clumps 
of diffuse emission (judged by eyes). 



Table 2. VNEI fitting results with the 90% confidence ranges and estimates of the gas density 
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fn e n H V/dl a (0.5-10 keV) 
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(10~ 2 ctss- 1 ) 




(10 22 cm~ 2 ) 
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(10 n cm- 3 s 


1 fin 57 ™ -3 ! no - 

I V CHI I I 1U 


-11 py.0-Q piyi 2 1 \ 
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3.88 ±0.08 


90.5/65 


2.8 ±0.1 


n fi7+ 002 

U-O'_o.04 


> 200 


9 cr+0.27 


0.84 




1 Q 
15 


2 


4.80 ±0.09 


87.0/79 


2.7 ±0.1 


^S+ - 06 

U.OO_ 05 


r i +4.6 
°- i -1.3 


3-84lJ;i 


1.31 




1 1 
1 1 


3 


2.49 ±0.07 


68.8/48 


3.4 ±0.2 


^fi+ 06 

U.OD_ og 


> 6.4 


o ^+2. 08 


1.05 




6.4 


4 b 


4.57 ±0.09 


108.9/79 


4 1+0-3 


62+ 04 
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9 42+3-10 

2 33 
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( [Mg/H]=1.12tg;Si 

9 7+ - 2 ^fi+0- 08 
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±S:Sg> [S/H]=0.80l°; 3 ) 








6 


3.30 ±0.08 


65.2/58 


5.ot2:5 


3 20+ 141 


0.97 




K 1 
O.l 


7 


3.31 ±0.08 
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o n+0.1 
J - u -0.2 


u - o4 -0.06 


> 5.5 


o Q7+1.65 
°- y '-0.98 


1.20 




D.o 


8 


3.08 ±0.07 


72.4/52 


2 9 +0 ' 2 

z - y -0.1 


0.63lg:8? 


> 3.4 


9 C9+O.6O 
Z - Oz -0.65 


0.95 






9 


2.23 ±0.06 


58.2/43 


2.9 ±0.2 


0.631HI 


> 2.7 


i-72lJ:iS 


0.62 




4.0 


10 


2.14 ±0.06 


44.2/44 




7Q+ 014 

u - ' y -0.10 


o 1 +3.6 
0.1_! 4 


1 ifi+o- 52 


0.40 




2.3 


11 


3.64 ±0.09 
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2.8 ±0.1 


0-5918:82 


> 3.8 
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J1 ' -0.77 


0.94 




1.9 

to 

2.0 ' 1 


12 


1.41 ±0.05 


49.8/32 
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0.58lS;S 
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3.1 ±0.1 
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Note. - 

a / denotes the filling factor of the hot gas. 

b Thawing metal abundances apparently improves fitting. 
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Table 3. Spectral fitting results for the NW unresolved source with the 90% confidence 

ranges 



regions net count rate 

(icr 3 ctss- 1 ) 


model 


xVd.c-.r. 


(10 22 cm- 2 ) 


photon index 


kT 
(keV) 


(0.5-10 keV) 
(10~ 14 ergscm^ 2 s _1 


5 2.60 ±0.20 
(1-10 keV) 


power 
bremsstrahlung 
blackbody 


4.08/4 
4.86/4 
6.95/4 


2 6+ 14 
6 +0 - 9 


o 9 + 1.3 


9 9+1.4 

0.75lg;i? 


9.1 
5.5 
3.3 
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5. Figure captions 



Fig. 1. — Panal a: Tri-color X-ray image of 3C 391. The X-ray intensity in the 0.3-1.5, 1.5- 
3.0, and 3.0-7.0 keV bands are color coded in red, green, and blue, respectively. The overlaid 
1.5 GHz radio contours are at 1.5, 3, 6, 9, 12, 18, 30, and 45xl0~ 3 Jy beam" 1 (from Moffett 
& Reynolds [1994]). Panel b: Tri-color X-ray + radio image of 3C391. The 1.5 GHz emission 
is plotted in red, the 0.3-3.0 keV emission is in green, and 3.0-7.0 keV is in blue. Panel c: The 
broad band (0.3-7.0 keV) X-ray image overlaid with 1.5 GHz radio contours (with the same 
levels described for the top panel). All the X-ray maps used here are exposure-corrected and 
are adaptively smoothed to achieve a S/N ratio of 3 (using the CIAO program CSMOOTH). 
The two cross labels in each panel denote the OH maser points (Frail et al. 1996). Panel d: 
Diffuse emission from SNR 3C 391 in the broad band 0.3-7.0 keV. The image is smoothed 
adaptively with a Gaussian kernel, which is adjusted to a count-to-noise ratio of 6. The 
color is scaled in the range from (1.19-331.85) x 10~ 2 ct s _1 arcmin -2 , logarithmically. The 
location of the sources removed from the data before the smoothing are marked by crosses. 
All the regions used for extracting spectra are indicated in blue, with cyan numerical labels. 
The dashed lines denote the border of the S3 chip. 

Fig. 2. — Chandra ACIS spectra of 12 small scale emission regions in 3C 391 fitted with the 
VNEI model. The region numbers, together with the factors multiplying the real fluxes, are 
labeled on the right side of the spectra. 

Fig. 3. — Chandra ACIS spectra of the medium-scale emission regions and entire remnant 
of 3C 391. The region numbers are labelled in each panel. 

Fig. 4. — Chandra ACIS spectrum of the NW unresolved source (region #5) fitted with the 
powerlaw model. 
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